Metagenomics is the study of all genomic content contained in given microbial communities. Metagenomic functional analysis aims to quantify protein families and reconstruct metabolic pathways from the metagenome. It plays a central role in understanding the interaction between the microbial community and its host or environment. De novo functional analysis, which allows the discovery of novel protein families, remains challenging for high-complexity communities. There are currently three main approaches for recovering novel genes or proteins: de novo nucleotide assembly, gene calling and peptide assembly. Unfortunately, their information dependency has been overlooked, and each has been formulated as an independent problem. In this work, we develop a sophisticated workflow called integrated Metagenomic Protein Predictor (iMPP), which leverages the information dependencies for better de novo functional analysis. iMPP contains three novel modules: a hybrid assembly graph generation module, a graph-based gene calling module, and a peptide assembly-based refinement module. iMPP significantly improved the existing gene calling sensitivity on unassembled metagenomic reads, achieving a 92–97% recall rate at a high precision level (>85%). iMPP further allowed for more sensitive and accurate peptide assembly, recovering more reference proteins and delivering more hypothetical protein sequences. The high performance of iMPP can provide a more comprehensive and unbiased view of the microbial communities under investigation. iMPP is freely available from https://github.com/Sirisha-t/iMPP.
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
more » « less
Abstract -
We performed various analyses on the taxonomic and functional features of the gut microbiome from NSCLC patients treated with immunotherapy to establish a model that may predict whether a patient will benefit from immunotherapy. We collected 65 published whole metagenome shotgun sequencing samples along with 14 samples from our previous study. We systematically studied the taxonomical characteristics of the dataset and used both the random forest (RF) and the multilayer perceptron (MLP) neural network models to predict patients with progression-free survival (PFS) above 6 months versus those below 3 months. Our results showed that the RF classifier achieved the highest F-score (85.2%) and the area under the receiver operating characteristic curve (AUC) (95%) using the protein families (Pfam) profile, and the MLP neural network classifier achieved a 99.9% F-score and 100% AUC using the same Pfam profile. When applying the model trained in the Pfam profile directly to predict the treatment response, we found that both trained RF and MLP classifiers significantly outperformed the stochastic predictor in F-score. Our results suggested that such a predictive model based on functional (e.g., Pfam) rather than taxonomic profile might be clinically useful to predict whether an NSCLC patient will benefit from immunotherapy, as both the F-score and AUC of functional profile outperform that of taxonomic profile. In addition, our model suggested that interactive biological processes such as methanogenesis, one-carbon, and amino acid metabolism might be important in regulating the immunotherapy response that warrants further investigation.more » « less
-
Noncoding RNAs (ncRNAs) play important regulatory and functional roles in microorganisms, such as regulation of gene expression, signaling, protein synthesis, and RNA processing. Hence, their classification and quantification are central tasks toward the understanding of the function of the microbial community. However, the majority of the current metagenomic sequencing technologies generate short reads, which may contain only a partial secondary structure that complicates ncRNA homology detection. Meanwhile, de novo assembly of the metagenomic sequencing data remains challenging for complex communities. To tackle these challenges, we developed a novel algorithm called DRAGoM (Detection of RNA using Assembly Graph from Metagenomic data). DRAGoM first constructs a hybrid graph by merging an assembly string graph and an assembly de Bruijn graph. Then, it classifies paths in the hybrid graph and their constituent readsinto differentncRNA families based on both sequence and structural homology. Our benchmark experiments show that DRAGoMcan improve the performance and robustness over traditional approaches on the classification and quantification of a wide class of ncRNA families.more » « less
-
more » « less
Abstract Motivation Accurate prediction and interpretation of ligand bioactivities are essential for virtual screening and drug discovery. Unfortunately, many important drug targets lack experimental data about the ligand bioactivities; this is particularly true for G protein-coupled receptors (GPCRs), which account for the targets of about a third of drugs currently on the market. Computational approaches with the potential of precise assessment of ligand bioactivities and determination of key substructural features which determine ligand bioactivities are needed to address this issue.
Results A new method, SED, was proposed to predict ligand bioactivities and to recognize key substructures associated with GPCRs through the coupling of screening for Lasso of long extended-connectivity fingerprints (ECFPs) with deep neural network training. The SED pipeline contains three successive steps: (i) representation of long ECFPs for ligand molecules, (ii) feature selection by screening for Lasso of ECFPs and (iii) bioactivity prediction through a deep neural network regression model. The method was examined on a set of 16 representative GPCRs that cover most subfamilies of human GPCRs, where each has 300–5000 ligand associations. The results show that SED achieves excellent performance in modelling ligand bioactivities, especially for those in the GPCR datasets without sufficient ligand associations, where SED improved the baseline predictors by 12% in correlation coefficient (r2) and 19% in root mean square error. Detail data analyses suggest that the major advantage of SED lies on its ability to detect substructures from long ECFPs which significantly improves the predictive performance.
Availability and implementation The source code and datasets of SED are freely available at https://zhanglab.ccmb.med.umich.edu/SED/.
Supplementary information Supplementary data are available at Bioinformatics online.
-
Janus nanoparticles (NPs) containing two chemically distinct materials in one system are of great significance for catalysis in terms of harnessing catalytic synergies that are not exist in either component. We herein present a novel synthetic method of two Janus-type MnOx-Ag and MnOx-AgI NPs for efficient photochemical water oxidation. The synthesis of Janus-type MnOx-AgI NPs is based on the oxidative nucleation and growth of Ag domains on MnO first and the subsequent iodization of Ag. A mild and non-disruptive iodization strategy is developed to yield Janus MnOx-AgI NPs, in which converting Ag to AgI domains with iodomethane (CH3I) is achieved through the partial iodization. Simultaneously, Mn2+ species in the primary MnO octahedrons are oxidized during the growth of Ag NPs, leading to the formation of amorphous p-type MnOx domains. Therefore, as-resultant Janus-type MnOx-AgI NPs combining two semiconductors into an integrated nanostructure can be used as an efficient photocatalyst for visible light-driven water oxidation. Janus MnOx-AgI NPs show an expected photocatalytic activity even in the absence of Ru(bpy)3Cl2 as an electron mediator. This intriguing synthesis may open up a new opportunity to develop asymmetric nanostructures of two semiconductors that will potentially be efficient photocatalysts for solar-driven water splitting.more » « less
